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1. Introduction 

The problem of random walks in a random environment (RWRE) has a long history and 
since the early results in the 1970s [1], a vast amount of informations have accumulated; 
for a recent review see Ref. [2j. The RWRE can be regarded as a toy model of 
disordered systems, for which exact results are available and which, due to its simple 
formulation, became a fundamental model in various fields such as transport processes 
or statistical mechanics of magnetic systems [3]. Most works concern the RWRE with 
nearest-neighbour jumps on the integers, for which a more or less complete picture 
is at our disposal. Beside rigorous results [U HI [5], this model was also studied by a 
strong disorder renormalization group (SDRG) method [6] which is closely related to 
that originally developed for disordered spin models [7]. This method, in which the 
small barriers of the energy landscape are successively eliminated, yields exact results 
for the asymptotical dynamics, among others the scaling of the typical displacement x 
of the walker with time t in the recurrent case: x ~ (\nt) 2 , in accordance with Sinai's 
theorem [I]. 

In higher dimensions, even on quasi-one-dimensional lattices or in case of non- 
nearest-neighbour jumps, the understanding of RWRE is at present far from complete. 
For the one-dimensional(lD) RWRE with bounded non-nearest-neighbour jumps, 
criteria for recurrence and transience are known [8] and for some special cases Sinai 
scaling was proven [9]. This model arises also in the context of disordered dynamical 
systems jXD] . For the RWRE on strips of finite width, which incorporates among others 
the former model and the persistent RWRE [11] , recurrence and transience criteria were 
obtained in Ref. [T2] . 

The aim of this paper is to propose an exact SDRG scheme for the RWRE on 
a strip. A necessary condition for the analytical tractability by the SDRG method is 
that the topology of the underlying lattice is invariant under the transformation, which 
generally does not hold apart from ID. As in our approach complete layers of lattice 
sites are decimated, the topology of the network of transitions is preserved. Contrary to 
the ID RWRE, the energy landscape does not exist in general, therefore we keep track of 
the transformation of jump rates in the same spirit as it was done for the closely related 
ID zero range process [13]. We shall show that in the fixed point, the transformation of 
relevant variables is identical to that of the ID RWRE with nearest-neighbour jumps, 
implying among others that Sinai scaling holds for strips of finite width in the recurrent 
case. 

The rest of the paper is organised as follows. In Section [2} the problem to be studied 
is defined in details. In Section [31 the renormalization group (RG) transformation is 
introduced and the RG equations are analysed in the recurrent well as in the 

zero-velocity transient phase. Finally, the results are discussed in Section HI 
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2. Formulation of the problem 

We consider a finite strip S — {1, . . . , L} x {1, . . . , m} of length L and width m, and 
call the set of sites (n, i) £ S with fixed n and i = 1, . . . , m the nth layer. We define 
on this lattice a continuous-time random walk by the following (nonnegative) transition 
rates for 1 < n < L: 



T{zt,z 2 ) 



P n (i,j) if z x = (n,i),z 2 = (n + l,j) 

Q n (i,j) if z\ = (n,i),z 2 = (n- 

if 2i = (n,i),z 2 = (n,j),i ^ j 
otherwise. 



Here and in the following, the formally appearing index (0, j) [{L + 1, j)] is meant to 
refer to site (L,j) i.e. the strip is periodic in the first coordinate. The m x m 

matrix P n (Q n ) contains the jump rates from the nth layer to the adjacent layer on 
the right(left), while the matrix R n with diagonal elements R n (i,i) := —J2j^iRn(hj) 
contains the intra-layer jump rates. Besides, we define the m x m matrix S n , which 
will be useful in later calculations by S n {i,j) := —R n (i,j), i ^ j, while the diagonal 
elements are fixed by 

(Pn + Qn~S n )l = 0, (1) 

where 1(0) is a column vector with all components 1(0). For the sequence of triples of 
matrixes, {(P n ,Q n , Rn)}, which defines the random environment, we impose at this 
point the only condition that it must be connected in the sense that every site is 
reachable from every other site through sequences of consecutive transitions with positive 
rates. The probability that the walker resides on site (n, i) in the stationary state is 
denoted by n n (i) and these are normalized as J2( n ,i) n n(i) = 1- Following Ref. [T2] , 
we introduce the row vectors ir n = {^n{i))i<i<m and for a fixed environment, write the 
system of linear equations that the stationary probabilities satisfy in the form: 

ir n S n = ix n -\P n -\ + ir n+1 Q n+ i, 1 < n < L. (2) 

Although, we started from a continuous-time random walk, the same equations can 
be written for a discrete-time jump process with transition probabilities obtained by 
rescaling the transition rates by max( n i ) S n (i,i). 



3. Renormalization group transformation 

The elementary step of the renormalization group method we apply is the elimination 
of the kth layer, such that the walker then jumps from the k — 1st layer directly to the 
k + 1st one with transition rates Pk-i{i,j) and from the k + 1st layer to the k — 1st 
one with rates Qk+i(i, j)- We choose the matrices Pk-i and Qk+i in such a way that 
the remaining L — 1 equations in ([2]) are fulfilled by the unchanged vectors ir n , n ^ k. 
Eliminating 77^ in Eq. ([2]), it turns out that also the matrices Sk-i and Sk+i must be 
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changed, and we have the following transformation rules: 




Sk+i — Sk+i — Qk+iS k Pk- 



(3) 
(4) 
(5) 
(6) 



All other matrices remain unchanged. The matrix S n has the following important 
property: 



which is meant to hold for the matrix elements. This can be proven as follows. We 
introduce the notation D m = det S n where the index m refers to the order of the matrix. 
The non-diagonal elements of S n are nonpositive, while S n (i, i) = J2j[Pn(h j)+Qn(h + 
J2j^iRn(hj) > for 1 < i < m since by assumption, the environment is connected. 
Regarding D m as a function of the variables e, := J2jS n (i,j) = J2j[Pn(h j)+Qn(h j)), i- e - 
D m = -D m (ei, . . . , e m ), it is clear that D m (0, . . . , 0) = and = D^-i where -D^_i is 
the determinant of the matrix S%> obtained from S n by deleting the ith row and column. 
Now, the relation D rn > can be shown by induction. Obviously, D 1 > 1. Assuming 
that -Dm-i = > ^ ^ or ^ — * — 171 an< ^ taking into account that connectedness implies 
J2i £i > 0, it follows that D m > 0. Thus det S n , as well as the diagonal elements of S' 1 
are positive. Using this result, the relations 5'~ 1 (i, j) > for i ^ j can be shown again 
by induction in a straightforward way. 

Relation ([7]) and Eq. ([5]) imply that ASk-i = Sk-i — Sfc-i = —Pk-iS^Qk < 0. In 
components: 



From these relations we obtain J2j AP k -i(i, j) < 0, where we have used AQ k -i = 0. 
Similarly, we obtain: AR k+ x(i, j) > 0, i ^ j and J2j AQk+i(i, j) < 0. Thus, the intra- 
layer transition rates are non- decreasing, while the sum of rates of inter- layer jumps 
starting from a given site is non-increasing under a renormalization step. 

Let us introduce the quantity Q n := I/US'" 1 !!, where the matrix norm || • || is defined 



as \\A\\ := maxj J2j From Eq. (jSJ), we have S k \ = S k \ + S k \P k -iS k 1 Q k S k l 



As relation (J7J) is valid also for the renormalized matrices, i.e. S k \, > 0, 

both terms on the right hand side are nonnegative, therefore ||*S r ^l_ 1 |j = + 
S^iPk-iS^QkS^^l > or, equivalently, £lk-i < fifc-i- By a similar calculation 

we obtain that Qk+i < ^fe+i- The RG procedure for finite L is defined as follows. 
The layer with the actually largest Q n is decimated, which results in a RWRE on a 
one layer shorter strip with effective rates given by Eqs. (JHHHJ) and the remaining ir n 
unchanged. This step is then iterated until a single layer is left. The variable defined 
by Q := max n n„, where n runs through the set of indices of non-decimated (or active) 
layers, decreases monotonously in the course of the procedure. For the special case 



(7) 



Ai2 A _ 1 (z,j)>0 (i^j) 



AS fc _i(i,i) < 0. 



(8) 
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m = 1 (ID), Qfc = Qfc(l, 1) + -Pfc(l, 1) and the transformation rules reduce to 




P fc _i(l,l)P fc (l,l) 

" g fc (i,i) + p fc (i,i)' 




Q fc+ i(i,i)g fc (i,i) 

" Q fc (l,l) + P fc (l,l)' 



(9) 



which have already been obtained in the context of the zero range process [T3] . 

The procedure described so far applies to any connected environment; as a trivial 
case even to the homogeneous environment. From now on we assume that the triples 
(P n , Qn, R n ) are independent, identically distributed random variables. We consider an 
infinite sequence of triples {(P n , Qn, R n )} and, in the usual continuum formulation [15] 
of the above RG procedure, we are interested in the asymptotic scaling of Q with the 
length scale £n that is given by the inverse of the number density cq of active layers: 



3.1. Recurrent case 

First, we focus on the case of transition rate distributions for which the random walk 
is recurrent in almost every environment. The question of recurrence is in general 
non-trivial for m > 1 (H [12]; nevertheless, a sufficient condition of recurrence is that 
the distribution of jump rates is invariant under the interchange of P n and Q n [14J. 
Furthermore, we do not deal with special environments which lead to normal diffusive 
behaviour (e.g. the case P n = Q n for all n). Instead, we consider less restricted 
situations: for instance, distributions where P n and Q n are independent. In this case, 
the above special environments form only a zero-measure set in the limit L — > oo. 

As a first step, we investigate the limits of transition rates when the density of 
active layers goes to zero. Consider a site (n, i) in an active layer in an arbitrary 
stadium of the RG procedure and assume that the initial matrix elements S n (i,j) 
were renormalized to some S n (i,j) < S n (i,j). Then we can write J2j^i R n (h j) < 
J2j^iRn(i,j)+T,j[Pn(i,j)+Qn(i,j)} = S n (i,i) < S n (i,i). Consequently, the intra-layer 
rates remain bounded throughout the RG procedure. Writing, e.g., eq. ([S]) in the form 
ASk-i = —Pk-iS^ Qk, we see that at least one of the sets of matrices {P n } and {Q n } 
must tend to zero as cn — > 0, otherwise the matrices S n would not remain bounded. 
Furthermore, it is clear that the assumption on recurrence requires that both {P n } and 
{Qn} must tend to zero if cq — > 0. This also implies that, in that limit, det S n — > 
and Q — ► Q* = 0. So, as the RG transformation progresses the inter-layer rates at the 
non-decimated layers are approaching zero without limits. 

For the study of various quantities close to the fixed point fl* = 0, it is expedient 
to define the following relation: / ~ g if lim^^Q f/g = 1. According to the above, 
we have Sk-i — Sk-i and similarly, for the matrix iS" 1 := S^/HS"" 1 !!, — <S A Tj L 1 

holds. One can easily show that the rows of 5" 1 are asymptotically identical, i.e. 
< ^n 1 ( i ')i) — ^n 1 (^)i) f° r 1 — hjfk < m, and the vectors formed from the rows 
tend to the stationary measure 7r„ of the isolated nth layer, i.e. 5~ 1 (z,j) ~ ^n(j) 



Cn = I/en. 
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for 1 < i, j < m, where ir n is the solution of the equation n n R n = which fulfils the 
condition Sj7f n (i) = 1. Although, the layers were not assumed to be connected within 
themselves initially, after many decimations they become almost surely connected due 
to the generated positive intra-layer transition rates when eliminating adjacent layers. 
If it is the case, the measure n n is unique. Introducing the matrices V n := S^ x P n 
and Q n := S~ x Q n) Eq. ([3]) can be written as Vk-i — Pk-i^k-i = Pk-iPk/^k with 
A k = S k \ - S k \. Using Eq. © we obtain that WS^Pk + Q fc )|| = 1. The rows of S k 1 
are asymptotically identical, therefore ||S'^ 1 Pfe|| + ||S'^ 1 Q/ C || — \\£>k (Pk + Qk)\\ = 1 
and Q k — W'PkW + IIQfell- Furthermore, A*. — * if Q — > 0, thus we obtain the 
asymptotical renormalization rule Vk-i — Pk-iPk/ (\\Pk\\ + 1| Qk\\), and we have a similar 
equation for Qk+i- Using that the rows of [jT-^H are asymptotically identical, we have 
H'Pfe-i'Pfcll — H'Pfc-ill • H'PfcH and obtain finally: 

nf> n ll^fc-ill " ll^ll ma H ||2fc+i|| ' IIQfcll / in >, 

iw+iiaii - l|gw "" m+m - ( ' 

We see that these equations have the same form as those of the ID RWRE in Eq. 
(J9j). The physical interpretation of these results is clear. If Q -C 1, the effective inter- 
layer rates are much smaller than the effective intra-layer rates, thus the walker in the 
renormalized environment spends very long time in a layer until it jumps to another 
one, so that its quasistationary distribution within the layer is given asymptotically by 
7r n . When the walker leaves the layer it does not "remember" at which site it entered 
the layer and irrespectively of this site, the effective jump rates to the adjacent layer 
to the right and left are \\V n \\ and ||Q n ||) respectively. Thus we may say that the 
model under study renormalizes asymptotically to a ID RWRE. In the course of the 
RG transformation, the normalization of the measure is obviously not conserved, i.e. 
S( n ,i) K n {i) < 1, where the prime denotes that the summation goes over the active sites. 
Nevertheless, on a finite strip, the walker spends most of the time in a small number of 
layers and the sum of ir n (i) over almost all sites goes to zero in the limit L — > oo, which 
is closely related to the Golosov localization [5]. At any stage of the RG transformation, 
the layer with the maximal Q n is decimated and Q n J2i 7i"n(*) can be interpreted, at least 
close to the fixed point, as the probability current from the nth layer to the neighbouring 
ones. This ensures that layers with smaller I^7r n (i), i.e. where the walker can be found 
with a smaller probability, are decimated typically earlier in the course of the SDRG 
procedure. Thus, fixing the length scale £ > 1 and renormalizing a finite strip of length 
L > £ to a strip of length L' = L/£, we expect that J2{ n ,i) ^rJS) ~^ ^(1) almost always if 
L — > oo. Now, if the correct normalization of 7r n (z) in the renormalized strip is restored 
by dividing by Z)( n ,t) n n(i), the probability current along the strip is modified only by an 
(9(1) factor. On the other hand, the current is invariant under the RG transformation, 
thus assuming that £ ^> 1, the RWRE on a strip of length L has the same current up 
to an 0(1) factor as an effective ID RWRE of length V ~ L. This implies that the 
current of the RWRE on a strip must asymptotically scale with the size as that of the 
ID RWRE. Consequently, the inverse of the current, which gives the mean time r that 
the walker needs to make a complete tour on the strip, must scale with L asymptotically 
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just as in one dimension: 

(lnr) 2 ~L. (11) 

Next, we have a closer look on the RG equations ( ITUl) and determine the scaling 
relation between Q and £n by pointing out the asymptotic equivalence to an already 
solved problem. In order to do this, we assume that the distributions of effective 
rates \\V\\ and ||Q|| broaden on logarithmic scale without limits as Q — ► 0. This 
property, which can be justified a posteriori, is characteristic of the so-called infinite 
randomness fixed points and ensures the asymptotical exactness of the procedure |15j . 
As a consequence, at the layer to be decimated, almost surely either H'Pfcll/HQfcH or 
II Qfc||/||Pfc|| tends to zero if Q — > 0. In the first case, Q ~ \\Vk\\ + ||Qfe|| — ||Qfc|| an d the 
decimation rules read 

Wh-iW ^ l|Pfc "l 'H Vk K \\Qk+i\\ ^ \\Qk+i\\, (12) 
II *ik\\ 

while in the second case Q ~ ||Pfc|| an d 

un-iii - \\Vk-ti \\Q k+ i\\ * llQk+ 'l' } Vhl] . (13) 

1 k\\ 



For the above transformation rules, it has been shown in Ref. [15] in the continuum 
limit that the distributions of \\V\\ and ||<2|| flow in the recurrent case (apart 
from some singular initial distributions) to the strongly attractive self-dual fixed 
point with identical distribution of \\V\\ and ||Q||: P*ijf) = e~ v Q(r]), where r\ = 
ln(Q/\\V\\)/ hx(p,o/Q), Qo is the initial value of Q and Q(x) is the Heaviside step function. 
Furthermore, the asymptotic scaling relation between £q and Q reads: 

^ ~ In 2 (fi /fi) . (14) 

Carrying out the RG transformation in a finite but long strip up to the last layer which 
is indexed by I, the magnitude of the current can be written as | J\ = \tti(Pi — Qi)\ ~ 
Si ^(OKII^II — IIQ/IDI ~ Ei 71 "^)^) where we used in the last step that for large L, 
\\Vi\\ and ||Q/|| differ typically by many orders of magnitude. Taking into account that 
is expected to remain finite for almost all environments in the limit L — > oo 
and substituting L for the length scale in Eq. ( Tl4l) we arrive again at Eq. ( fTTT) . From 
this scaling relation we conclude that the typical displacement of the first coordinate x 
of the walker on an infinite strip scales with the time in the recurrent case as x ~ (lnt) 2 
for almost all environments. 



3.2. Sub-linear transient phase 

Now, we consider the case when the environment is still an independent, identically 
distributed sequence but the random walk is transient. It is known for the ID 
RWRE that if < /^i < 1, where /ii is the unique positive root of the equation 
[Q(l, 1)/P(1, 1)] M1 = 1 and the over-bar denotes averaging over the distributions of 
(5(1,1) and P(l,l), the displacement grows sub-linearly as x ~ t Ml [TJ [16]. In the 
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analogous zero-velocity transient phase of the RWRE on a strip, the matrices P n and 
Q n must still renormalize to zero, and the asymptotical transformation rules are given 
by Eqs. ( fl2lfT3l) . The analysis of these RG equations in the continuum limit has been 
carried out in Ref. [17] and has yielded the asymptotical result: £q ~ (f2/fi ) _M . We 
thus conclude that the displacement grows as x ~ t^ also for the RWRE on a strip in this 
phase. For the ID RWRE, fj, — y^, which is due to the fact that the energy landscape 
defined by U n+ i — U n = \n[Q n+ i(l,l)/ P n (l,l)] carries the full information on fii and 
even the approximative rules in Eqs. ( fl2iri3i) leave the energy difference between active 
sites invariant (cf. the method in Ref. [6]). For m > 1, Eqs. (I12H13P are valid only 
asymptotically and the problem how the exponent /i is related to the initial distribution 
of jump rates is out of the scope of this approach. 

4. Discussion 

We have presented in this work an SDRG scheme for the RWRE on quasi-one- 
dimensional lattices, which incorporates also the RWRE with bounded non-nearest 
neighbour jumps. We have made use of that by eliminating appropriately chosen 
groups of lattice sites, the topology of the network of transitions remains invariant. 
We mention that there are special sub-networks of transitions with positive rates which 
are invariant under the transformation: As can be seen from Eqs. ([3]|I]), if the ith row 
or column of P n or Q n is zero for all n, then this remains valid also after an RG step. 
An example for m = 2 is the process with the only positive inter- layer rates P n (l, 1) 
and Q n (2, 2), which can be interpreted as a ID persistent RWRE. We have shown that 
the model renormalizes to an effective ID RWRE and concluded that, although, the 
finite-size corrections are strong (see Ref. pH]), Sinai scaling is valid asymptotically in 
the recurrent case, while in the sub-linear transient regime the displacement grows as 
x ~ t M . Although, the method is not appropriate for establishing an analytical relation 
between the non-universal exponent fi and the initial distribution of jump rates, the 
numerical implementation of the exact RG scheme provides a more efficient tool for the 
estimation of /i than the direct solution of Eqs. (T5]). 

When this work was finalized, a preprint by Bolthausen and Goldsheid appeared, 
in which similar results are obtained in the recurrent case in a different way [IS] . 
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